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The double pulsar J0737-3039A/B is a unique system with which to test grav- 
itational theories in the strong-field regime. However, the accuracy of such 
tests will be limited by knowledge of the distance and relative motion of the 
system. Here we present very long baseline interferometry observations which 
reveal that the distance to PSR J0737-3039A/B is llSOligo PC» more than dou- 
ble previous estimates, and confirm its low transverse velocity (~9 km s^). 
Combined with a decade of pulsar timing, these results will allow tests of grav- 
itational radiation emission theories at the 0.01 % level, putting stringent con- 
straints on theories which predict dipolar gravitational radiation. They also 
allow insight into the system's formation and the source of its high-energy 
emission. 

The double pulsar system J0737-3039A/B (I7]|2j) is one of eight known double neutron star 
(DNS) systems, and the only system in which both neutron stars are visible as pulsars. The 
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'A' pulsar, which has been spun-up by accretion ("recycled") to a period P = 22.7 ms, was 
discovered first; the non-recycled 'B' pulsar (P = 2.77 s) was found during follow-up timing 
observations. PSR J0737-3039A/B is the most relativistic known DNS system. It has an orbital 
period of 2.5 hours and a coalescence time (due to orbital energy loss to gravitational radiation) 
of 85 Myr. Compared to most DNS systems, it has a low eccentricity (0.08) and is thought 
to possess a low transverse velocity [10 km s~^; ©], which is difficult to explain in standard 
models of pulsar formation. Pulsar timing of PSR J0737-3039A/B produces results consistent 
with the theory of general relativity (GR) at the 0.05% level ©. 

The distance to PSR J0737-3039A/B is estimated using the pulsar dispersion measure (DM) 
and models of the ionized component of the Galaxy (I?]©. However, distances derived in this 
manner have been shown to be in error by a factor of two or more for individual systems [e.g. 
PSR B0656+14; I©]. A more accurate distance to the PSR J0737-3039A/B system is needed 
in order to determine the contribution of kinematic effects to pulsar timing, whose uncertainty 
would otherwise limit the precision GR tests could achieve; to establish the source of the x-ray 
emission from the system ([2 HI); and to refine the estimated radio luminosity. 

We used the Australian Long Baseline Array (LBA) to obtain a model-independent dis- 
tance to PSR J0737-3039A/B, through the direct measurement of its annual geometric parallax. 
We have made seven very long baseline interferometry (VLBI) observations of PSR J0737- 
3039A/B over an 18 month period between August 2006 and February 2008, which allowed us 
to fit the system's position, proper motion and parallax [for more details see Table [T] and 
Fig.m all errors are la unless otherwise stated]. 

The distance to PSR J0737-3039A/B is llSO^igo P'^' which is inconsistent with previous 
DM-based estimates of 570 pc (4) and 480 pc (5). Our measurement is more significant than 
the distance of SSSiill pc obtained through a marginal timing parallax detection (J), which 
differed by approximately la from our result. The revised value of electron density along the 
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line of sight to PSR J0737-3039A/B is 0.043 cm"^ suggesting that the influence of the Gum 
nebula (a large, ionized hydrogen region between the Earth and PSR J0737-3039A/B) along 
this line of sight is less than originally thought. The proper motion of PSR J0737-3039A/B is 
4.37 ± 0.65 mas yr^^, consistent with a previous timing measurement of 4.2 ±0.6 mas yr~^ 

We used the revised distance to PSR J0737-3039A/B to estimate the precision attainable 
for testing predictions of gravitational wave emission from the system. In order to compare the 
rate of change of orbital period (Pb) due to the loss of energy to gravitational radiation with 
the GR prediction, the observed Pb must be measured as accurately as possible [the current 
significance is 70 cr (3)], and contributing factors to Pb other than GR must be estimated and 
subtracted as accurately as possible. For PSR J0737-3039A/B, the major contributing factors 
are differential Galactic rotation (Pb°*), acceleration towards the plane of the Galaxy (Pb), and 
the apparent acceleration caused by the transverse motion of the system [the Shklovskii effect 
Pb^"^; (|7Q|)1. to which we will refer collectively as galactic and kinematic contributions (Pb'^). 
For the newly-calculated distance of 1 150 pc, the magnitude of these effects is calculated using 
Equations 2.12 and 2.28 from (|77)) as: 

Pb cRo \ sm I + (j^ J 

— = sm6 (2) 

Pb c 

^ = ^ (3) 

Pb c 

where /, h and z are galactic longitude, latitude and height respectively, d and fi are pulsar 
distance and proper motion, Rq and vq are the Galactic radius and speed of the Solar System 
[taken to be 7.5 kpc and 195 km respectively ([721)1. is the vertical gravitational potential 
of the Galaxy [taken from ^ as 0.45 km^ s"^ pc^ at the height of PSR J0737-3039A/B] and 
c is the speed of light. The dominant uncertainties are in d (16%), fi (15%), and Rq, vq and K^, 
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whose errors are estimated at ~10%. 

Equations©-® give Pb'°V^b = (-4.3±0.7) x IQ-^o §-1, P^7Pb = (3.8±0.8) x IQ-^i s-\ 
and P^^^/Pb = (5.3 ± 1.8) x lO^^o s^^. Combining these terms gives P^^/Pb = (1.5 ± 2.1) x 
10~^° s^^, and multiplying by the observed orbital period P^ = 8834.5 s (3) yields the net effect 
of these terms on the observed orbital period derivative: P^^ = (1.3 ± 1.8) x 10^^^. 

These contributions to the orbital period derivative are four orders of magnitude below the 
GR contribution, and two orders of magnitude below the current measurement error [P{^^^ = 
( — 1.252 ± 0.017) X 10~^^; Qi]. Thus, with the current accuracy in the measurement of distance 
and transverse velocity, GR tests can be made to the 0.01% level with PSR J0737-3039A/B 
using Pb. However, about ten years of further precision timing will be required to reach this 
point. Measuring Pb at this level will place stringent requirements on the class of gravitational 
theories which predict significant amounts of dipolar gravitational radiation, exceeding the best 
Solar System tests ©. Measuring the moment of inertia of pulsar A, however, would require 
another order of magnitude improvement in the measurement precision of Pb Q. In the near 
future, additional VLBI and/or timing measurements can be expected to reduce the error in both 
the distance and velocity of PSR J0737-3039A/B below 10%; however, even with negligible 
error in these parameters, the existing accuracy of measurements of Rq, Vq and would limit 
the accuracy of Pb measurements in this system to 0.004%. To attain the 10~^ precision neces- 
sary to measure the neutron star moment of inertia, the constants Rq and fo must be measured 
to a precision approaching 1%. 

X-ray observations of PSR J0737-3039A/B show that most of the x-ray emission from 
the system is modulated at the spin period of the A pulsar (01241 ID 1221), but there has been 
considerable debate over where and how the x-rays are generated. Normally, pulsar x-ray 
emission is thought to have a magneto spheric or thermal origin (|76l) . However, the small binary 
separation and interaction of the pulsar wind of the A pulsar with the magnetosphere of the B 
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pulsar (1771) provides alternate x-ray generation mechanisms. 

The spin-down luminosity of pulsar A is over three orders of magnitude greater than that 
of pulsar B. Thus, pulsar A is likely to dominate any magneto spheric x-ray emission from the 
system. Whilst a wholly magneto spheric (power-law) origin for the observed x-ray emission 
is plausible, the neutral hydrogen column density [A^h ~ 1.5 x 10^^ cm~^; @] implied by the 
model is higher than expected from the pulsar's previously assumed location in the Gum nebula. 
The Gum nebula is believed to be approximately 500 pc distant, with a depth of several hundred 
pc ©. 

However, the value of Ah is consistent with the usual average of 10 neutral hydrogen atoms 
for every free electron along the line of sight. Our revised distance estimate places PSR J0737- 
3039A/B beyond the Gum nebula, implying that the measured value for A^h is not discrepant. 
It also increases the estimated x-ray luminosity by a factor of 5, but the revised value for a 
power-law fit (1.2 x 10^^ erg s^^) remains consistent with known relations between pulsar spin- 
down luminosity and x-ray luminosity OM. Hence, our result supports a power-law model of 
magneto spheric origin (from pulsar A) for the bulk of the x-ray emission from PSR J0737- 
3039 A/B. 

The discovery of PSR J0737-3039A led to a marked upward revision in the estimated Galac- 
tic merger rate of DNS systems (|79l) . although uncertainty over the characteristics of recycled 
pulsars means that the true value of the merger rate remains poorly constrained. Specifically, 
the distribution of recycled pulsar luminosities is generally extrapolated from the entire pulsar 
population (|79l) even though it (along with the distributions of pulse shape and beaming frac- 
tion) appears to differ from the distribution for slower pulsars (|2Q|) . Our revised distance shows 
that the radio luminosity of PSR J0737-3039A is a factor of five greater than previously as- 
sumed. If this revision were to markedly influence the recycled pulsar luminosity function, then 
the assumed space density of DNS systems would be reduced, with a corresponding impact on 
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DNS merger rates estimations. 

Finally, we used the measured transverse velocity for PSR J0737-3039A/B (24lgkm s~^) 
to constrain models of the formation of the system. After subtracting estimates of the peculiar 
motion of the Solar System and Galactic rotation (21), we measure a transverse velocity in 
local standard of rest of 9^3 km s^^. This is comparable to the unadjusted value of 10 km s^^ 
presented in (3 ), and is within the range of transverse velocities expected for the massive stars 
which are DNS progenitors [~20 km ; (|22))] . Because the transverse velocity of PSR J0737- 
3039A/B is so low, if the system received a large velocity kick at birth, it must have a large 
radial velocity. However, there are no observational methods available to determine the radial 
velocity in a DNS system. 

Because of the accurate measurement of its Shapiro delay, PSR J0737-3039A/B is known 
to lie edge-on (3). If the only kick it received was provided by the loss of binding energy 
during the supernova explosion, the resultant 3D space velocity should be on the order of ~50 
km s~^, estimated from the system's observed eccentricity and orbital velocity (|25l) . This space 
velocity would be constrained to the plane of the orbit. From simple geometry, the probability 
of observing a transverse velocity less than 10 km s^^ is about one in eight, which is small, 
but not unreasonable. Conversely, if the double pulsar had received a large kick (|2?|) . the odds 
of observing such a low transverse velocity become increasingly remote. Not only would the 
radial velocity have to be increasingly large, but the inclination angle of the system must not 
be altered by the kick. Hence, our transverse velocity results reinforce those of © and are 
consistent with the interpretation of (|25l) . who argue for almost no mass loss and kick in the 
case of PSR J0737-3039A/B. 

The implication of low kick velocities in PSR J0737-3039A/B-like systems offers a pos- 
sible, albeit speculative, explanation to the formation of PSR J 1903+0327, a heavy, highly 
recycled millisecond pulsar (mass 1.8 solar masses, period 2.15 ms) with a main sequence com- 
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panion of one solar mass (l26l) . The orbit of such a pulsar should have been highly circularized 
during the mass-transfer phase (|27l) . However, PSR J1903+0327 possesses an intermediate 
orbital eccentricity (e = 0.44). 

A formation mechanism for PSR J 1903+0327 has been suggested in which a triple system 
experiences a white dwarf-neutron star coalescence (28\). However, a coalescing DNS system 
such as PSR J0737-3039A/B could also create a PSR J1903+0327-like pulsar. Thus, given the 
low velocity of PSR J0737-3039A/B, an alternative formation mechanism for PSR J 1903+0327 
involves a triple system containing a close DNS binary and a main sequence star. 
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Table 1 . Fitted VLB! results for PSR J0737-3039A/B 



Parameter Value 

Right Ascension a (RA; J2000)^ 
Declination S (Dec; J2000)^ 
Proper motion in RA (/^ q, X cos 5 ; mas yr ~ "'" ) 
Proper motion in Declination (fi^ ; mas yr~ ) 
Parallax (mas) 
Distance (pc) 

Transverse velocity (km s ~ "'" ) 
Transverse velocity in LSR^(km s~ ^) 
Reference epoch (MJD) 



As discussed in (9], the actual error in the pulsar position is dominated by the 
alignment of the barycentric reference frame used for pulsar timing and the quasi- 
inertial frame used for VLBI, and is approximately an order of magnitude greater 
than the formal fit error shown here 

^ Local Standard of Rest 
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Fig. 1: Motion of PSR J0737-3039A/B plotted against time. (Top) Offset in right ascension; 
(Bottom) Offset in declination. The best fit (dashed line) is overlaid on the measured positions. 
The reduced chi-squared of the fit was 0.79, implying that the measurement errors (and thus 
the errors on fitted parameters) might be overestimated and hence conservative. 
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Supporting Online Material 



Materials and Methods 

Details of the observational, correlation, and post-processing strategies used are covered in 
detail in (SI) and are summarized here. Seven observational epochs spanned an 18 month pe- 
riod from August 2006 to February 2008. Approximately six hours of on-source time (divided 
equally between the phase reference calibrator B0736-303 and PSR J0737-3039A/B) were 
available per epoch, and all available antennas from the Australian Long Baseline Array (LBA) 
were utilized. The LBA is a network of six radio telescopes which provides the best means for 
VLBl study of southern radio sources. Of these, five telescopes are equipped to make observa- 
tions at 1650 MHz, which was the observing frequency utilized. These were the three ATNF 
antennas (ATCA, Parkes, Mopra), the University of Tasmania antenna near Hobart and, when 
available, the 70m NASA Deep Space Network (DSN) antenna at the Tidbinbilla station. In 
each session, the on-source time was spread over the ten hour period that PSR J0737-3039A/B 
was visible to the array in order to obtain good uv coverage. Four 16 MHz bands were used, 
quantized at two bit precision to give a total data rate of 256 Mbps. 

The phase reference calibrator B0736-303 is a bright and predominantly compact radio 
source separated by ~20 arcminutes from PSR J0737-3039A/B. Since the position of B0736- 
303 was not known precisely in the International Celestial Reference Frame (ICRF) at the outset 
of this project, its position was determined by aligning the obtained position for PSR J0737- 
3039A/B with the value known from timing (S2), accounting for the pulsar proper motion. This 
is possible due to the alignment at the milliarcsecond level of the ICRF with the barycentric 
reference frame used for pulsar timing. The final positional accuracy for B0736-303 was esti- 
mated to be better than 5 milliarcseconds, which is acceptable for the short calibrator throw (20 
arcminutes) used here. This use of the timing position of PSR J0737-3039A/B to obtain the 
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calibrator position means that the VLBI observations do not contribute to an improvement in the 
reference position for PSR J0737-3039A/B, as the error in the ultimately obtained position is 
dominated by the uncertainty in alignment between the ICRF and barycentric reference frames 
[1-2 milliarcseconds; (S3)~\. A six minute target/calibrator cycle was employed, with observing 
time divided equally between the target and calibrator. 

The data were correlated, using matched filtering (gating) on pulse profiles, with the DiFX 
software correlator (S4). Pulse profile filtering improves the signal to noise when observing 
a periodic signal by weighting the correlator output according to the varying pulsar flux. For 
PSR J0737-3039A/B, the pulse filter was set based on the signal from pulsar A only, since this 
dominates the emission in the system. The pulse filtering yields a factor of 2.5 improvement in 
the observed signal to noise ratio. The most accurate published pulsar ephemeris (S2) was used 
to set the pulse filter. The visibility data for each epoch was weighted by baseline sensitivity 
and imaged using natural weighting. 

The data were fitted using an iterative approach designed to estimate and account for sys- 
tematic errors, described in detail in ( SI ). This approach generated a position and position error 
for each 16 MHz band at each epoch, and the single-band positions were used to obtain a single 
weighted centroid position for each epoch. Estimates of systematic error at each epoch were 
made by inspecting the reduced chi-squared of the fit to position centroid, and these estimates 
were added in quadrature to the centroid position error. The final seven position measurements 
were then fit for position, proper motion and parallax using the iterative minimization code de- 
scribed in (SI). The reduced chi-squared of the final fit was 0.79, indicating that the quoted 
errors for the fitted values may be slightly conservative. 
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